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Abstract 

A neutrino produced in a particle decay maintains an unusual wave nature that reveals a diffraction phenomenon in a 
large spatial region. Diffraction gives a new component to a neutrino flux in this region and modifies determinations 
of physical quantities in neutrino experiments. A total cross section of neutrino nucleon scattering at high energies 
is studied in this paper. A cross section is proportional to an energy in the quark-parton model and is modified with 
the energy-dependent neutrino flux induced by the diffraction component. A total number of events devided by the 
energy is not constant but has a slight energy dependence. Theoretical values agree well with the recent experiments. 
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1. Neutrino diffraction and absolute neutrino flux 

In high energy neutrino experiments, neutrinos from 
pion decays are used. Neutrino flux is not directly mea- 
surable and is estimated from a pion's flux and decay 
probabilities. Now precision measurements are becom- 
ing possible in neutrino experiments and a precise value 
of the flux is required. So it is important to verify a stan- 
dard assumption and to find a precise value. We point 
out that a correction due to a quantum mechanical inter- 
ference of a single neutrino is necessary in the flux in 
near detector regions. 

Decay amplitudes and probabilities are obtained nor- 
mally with a standard method of S-matrix using plane 
waves. An asymptotic boundary condition HI 12 is as- 
sumed and an average decay rate is found and used for 
estimations of the flux. This method gives precise val- 
ues of fluxes of final states in the asymptotic space- 
time region where a particle has no correlation and is 
described with a plane wave. It has never been clari- 
fied, however, if particles in the final state are observed 
in an asymptotic region. It would be reasonable to as- 
sume that a boundary between an asymptotic and non- 
asymptotic regions is of microscopic size, since tran- 
sitions of weak decays occur in short distance regions. 
However we find that is not the case. An asymptotic 
region for a neutrino in a particle decay is located in 
a macroscopically distant area from an initial particle's 
position. Waves produced at different positions accu- 
mulate at a detector and make experiments in the inside 
of this area unusual. Superposition of waves causes an 



interference effect in a form of a diffraction in a paral- 
lel direction to the momentum, which makes the plain 
waves non-asymptotic. Various relations satisfied for 
the asymptotic states, such as a unitarity, an energy mo- 
mentum conservation, and others are modified with the 
diffraction term. 

To study the interference phenomenon of particles 
produced in a decay, it is convenient to study a coher- 
ence length from a behavior of a correlation function 
defined with a wave function. Let i[f(x,y) be a wave 
function of a final state composed of one particular par- 
ticle at a space-time position x and others, y. An am- 
plitude of observing this particle at a momentum p and 
a position X using a detector is defined with the wave 
function at one space-time position x and a wave packet 
w(p, x) in the form 



m{X, p, x, y) = w(p, x - X)\]/{x, y), 



(1) 



where w(p, x - X) is a wave packet with which this par- 
ticle interacts and a Gaussian form with a size o~ and 
velocity v = p/E, E — y/? 2 + m 2 is used for the sake of 
simplicity, 



w(p, x - X) = Noe 



-ip-(x-X) 



(2) 



where p ■ x is scalar product of a four dimensional co- 
ordinate x and momentum p and No is a normalization 
factor. 

Since the wave packet vanishes at \x\ — > oo, it sat- 
isfies the asymptotic boundary condition. So the wave 
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packet is approapriate to study a finite size correction. 
An integral of the product x V*(X,p,x\,yYV(X,p,X2,y) 
is the probability C(X, p), 



C(X, p) = j d\ l dyd A x 2 x V{X, p; , x u y^X, p; x 2 , y), 

where dy stands for the volume element of state y. 
C(X, p) is written as 

J d 4 x\d 4 X2W*(p, x\ - X)w(p, X2 - X)A(xi,X2), 

using the correlation function defined with 

A(xi,jc 2 ) = j dy4i(x\,y)*\p(x2,y). (3) 

The correlation function generally has the singular func- 
tion 



6(A)e(x° 1 - x° 2 ),A = (jci -x 2 f, 



(4) 



or regular functions e" nA , e~ mA where a finite m is deter- 
mined from dynamics of the system. If one momentum 
that is conjugate to x\ - X2 is integrated and the singular 
function 6(A) is derived, this gives the form 



C(X,p) = C m e 'o, |f| > l , 

where Iq is given in the form 

2\p\tw 
k = 5— , 



(5) 



(6) 



in a high energy region, \p\ » m. Other oscillating or 
decreasing terms give constants to C(X). So Iq stands 
for a coherence length of the wave function. The coher- 
ence lengths of a pion and an electron of an energy of 1 
[GeV] are 



[GeV- 1 ] = 2x KT 14 [m], (7) 
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z ,/™„ = __ [GeV !] ^ 10 10 [m], (8) 

and those of other hadrons are shorter than that of a 
pion. So charged leptons and hadrons have short coher- 
ence lengths of microscopic sizes in an ordinary high 
energy region. Neutrinos are exceptional and lighter 
than an electron by 10 6 or more and have coherence 
lengths, 



jmutrtno w jq2 - IQ 3 [m], 



(9) 



which is a macroscopic size. From Eqs. (O and (0, the 
non-asymptotic region is narrow in charged leptons and 
very wide in neutrinos. 



In outside of the coherence length \X\ » Iq, the wave- 
like correlation vanishes. Hence the state behaves like a 
particle-like object which has no spatial correlation and 
satisfies the asymptotic condition. In inside, \X\ < Iq, 
on the other hand, the probability varies with |a] and 
a wave-like correlation remains. So the physical state 
behaves like an unusual wave that is distinct from a 
particle-like object . The physical state is not treated as 
a real particle and does not satisfy the asymptotic con- 
dition. So the coherence length is used as the boundary 
between the asymptotic and non-asymptotic regions. 

Transition probabilities in the asymptotic region are 
studied using an ordinary S-matrix and a decay ampli- 
tude of plane waves and hold relations of a standard 
S-matrix. But those in the non-asymptotic region are 
peculiar and different from those of particle natures of 
a standard S-matrix. They violate some relations and 
are not studied with a standard S-matrix of plane waves. 
The physics in both regions can be studied with a time 
dependent Schroedinger equation or a scattering am- 
plitude of wave packets, which satisfy the asymptotic 
conditions and have manifest position dependence. We 
study the physics in the non-asymptotic region using the 
position dependent amplitudes. 

Charged leptons are massive and their measurements 
are made in their asymptotic regions. So physical quan- 
tities are computed with the ordinary S-matrix and agree 
with values obtained in experiments. Weak decays of 
mesons are understood well with V - A weak interac- 
tions through measurements of the charged leptons. For 
neutrinos the situation is different. Since the asymp- 
totic region is located in a distant area, neutrinos may 
be observed in non-asymptotic regions. Then the neu- 
trino reveals wave-like phenomena such as interference 
or diffraction and has an unusual position dependence. 
A diffraction term, in fact, emerges and gives a finite 
contribution to the neutrino flux. A new term which 
does not satisfy varios relations of the standard S-matrix 
is added to an amplitude in a short distance region. The 
total neutrino's flux thus obtained deviates from that ob- 
tained using the naive decay probability and has an ex- 
cess. Implications of the diffraction term to the total 
cross sections [3, 4,. 3] are studied. 



2. Position dependent probability 

A diffraction term 10, 01 is derived from a position- 
dependent amplitude of a pion decay process. The am- 
plitude is determined in the first order of the V -A weak 
Hamiltonian H H . in the form, T — J d 4 x{l,v\H w (x)\n), 
here a pion is prepared at a time T„, and a neutrino 
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is measured with a wave packet of a space-time posi- 
tion (T v ,Xy) and a charged lepton is un-measured. The 
neutrino wave packet (a |9l 1 1011 expresses a target nu- 
cleon which the neutrino interacts with. Hence the wave 
packet is well localized in the coordinate variable. To 
represent this property, the momentum variables must 
cover whole momentum region and the Gaussian form 
is used in this work. A set of wave packets of a cr v 
with a continius spectrum of momentum and coordi- 
nate is complete ja]. Hence neutrinos are described 
with wave packets of central values of the momentum 
and coordinate and a width ITll[H[l3i[l4l[l5i[l6[ [l7l. 
Other particles are described with the plane waves. 
They are expressed in the form \n) - \p„, T„), \l, v) = 
\l,Pi; v, p Y , X V ,T V ). The amplitude T is written with the 
hadronic V - A current and Dirac spinors in the form 

T = j d 4 xdk v N 1 (0\J^_ A (x)\7T) 

x Hpdy.d - r5 )v(^ y ) e ^-+^-(-Xv)-¥(?v-A) 2 ) (10) 

where Ni - ig((r v /jf) J (mim v / 'EtE v )- , and the time t is 
integrated in the region T^ < t. cr v is a neutrino wave 
packet size and is estimated using a nucleus size. The 
Gaussian form of the wave packet is used for a sake of 
simplicity but the most important result is unchanged in 
general wave packets as was verified in Jg] . 

If the coordinate x is integrated in Eq. (II Oi l, the 
delta function (2n) i 6{p n — pi — k v ) emerges. So k v 
is integrated easily and the Gaussian part becomes 
e ~it(?n-pi-Pvf +iE(p n -pi-p v )(/' '-T v )-i(p„-pi-p v yx v _ rr; ne eX p . 

nent has two stationary momenta, pi ~ p n - p v and 
- T v ) - % ~ OJ = a J&zjjzM H. Hence the 
probability gets a contribution from a lepton around the 
former momentum, which is called a normal term, and 
another contribution from it of broad spectrum around 
the latter momentum, which is called a diffraction term. 

The transition probability to this final state is finite 
and an order of integrations are interchangeable. So 
the probability is written, after the spin summations are 
made, with a correlation function and the neutrino wave 
function in the form 



f dpi y 

EyJ 



51,52 

d A x x d A x 2 e-^ ^-^ )2 A nl (6x)e mx \ (11) 
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where N 2 = g l (4n/o- v )i V~\ V is a normalization vol- 
ume for the initial pion, xf — X v + v v (f,- - T v ), 5x = 



x\ — X2, <p(Sx) — p v -5x and 

K,i(5x) =7^-r \ -^r(2p„ ■ p v pn ■ pi - m\pi ■ p v ) 
(2nY J E(pi) 

Xe -i(p»-Pi>Sx_ (12) 

In Eq. (fT2l i. the momentum, pi, is integrated in whole 
positive energy region. 

3. Light-cone singularity 

A K i(6x) is composed of a light-cone singularity 
|@, [3 and regular terms. The former is generated 
from those plane waves that have a same phase and 
are added constructively. In order to extract a leading 
singular term in the variable, 6x, we write the inte- 
gral in a four dimensional form with a new variable 
q = pi — p„ that is conjugate to 5x. Then A„j(6x) is 
decomposed into integrals in < q° and -p\ < q° < 0. 
An integral from < q° is written in the form, 
[2{p n ■ p v )p n ■ [p n - igfj) - ml (p n - / Jj) • p v ] I u 
where 



4tt 4 



1 



q 2 + 2p„-q + m 2 — ie 



Jq-Sx 



and m = m n - mf. The integrand of I\ is expanded in 
p n -q and the integration leads the light-cone singularity 
01811 . 5{5x 2 ), and less singular and regular terms which 
are described with Bessel functions. An integral from 
the region — p\ < q° < 0, h, is written with the mo- 
mentum q — q + p„ and has no singularity. Thus the 
correlation function, A K j(Sx), is written in the form 



A nJ (6x) = 2i < 2{p K ■ p v )p„ \Pn-i 



d5x 



Di„ I -i^r) (^^) + f short ) + /.: 



(13) 



where A = 6x 2 , D*(-i^) = Zi (l/l\)(2 P A-i£) ■£?)', 
f shon = e(-A) [Ni(&- i e(5t) MZ)}- Wi e(A) 
K!(g), £ = m VI, NuJu and K\ are Bessel functions. 
f short has a singularity of the form 1 1 A around A — and 
decrease as e~ m ^ or oscillates as e ' nl VUI a t a large \A\. 

The series in Eq. ( ITJb converges when 2p n -p v < 
and this expression is valid in this kinematical region. 

The singular terms, 6(A) and others, in Eq. ( fT3l de- 
crease slowly with the distance and give a correlation 
effect in a wide area. These terms are derived from 
the integration in E n < Ei. Because this region is out- 
side of kinematical region of satisfying the energy and 
momentum conservation p n — p v + p\, its contribution 
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disappears in the amplitude defined at T = oo. Con- 
versely the singular terms are included only in the phys- 
ical quantities observed at a finite time interval. So they 
are not derived from the standard calculations of plane 
waves with the asymptotic conditions. The latter term, 
h, on the other hand, comes from the integration region 
Ei < E n , which is in the kinematical region of satisfying 
the energy and momentum conservation. So this deter- 
mines the quantities at T = oo. This term oscillates or 
decreases fast in A with a time scale determined with or- 
dinary microscopic quantities and becomes microscop- 
ically short, as most other cases. We will see that the 
light-cone singularity is combined with the small neu- 
trino mass and gives a finite distance correction of an ex- 
ceptional scale to the position dependent neutrino flux. 

4. Integration of space-time coordinates 

Next, Eq. (TTJt is substituted to Eq. (fTTT i and the coor- 
dinates x\ and X2 are integrated. The most singular term, 
Js(A), is from ^r~8(A) in A„j, which has no scale, and is 
rewritten using the center coordinate X M - (x^ + x%)/2 
and the relative coordinate r - x\- X2- The center coor- 
dinate X is integrated easily and Jg{x) becomes an inte- 
gral of the relative coordinates (rj, r{). Finally we have 



h(X) = C 5 a)e(dt)\dt\ l e ' ' 



(14) 



where Cg{x> = (ovr)icr v /2 and (f> c (5t) = w v 5t,u) Y = 
rriy/2E v . The phase 4>(6x) of Eq. (fTTT i becomes the small 
phase (f>c(5t) of Eq. (TBI at the light cone due to a can- 
cellation between the time and space components. The 
angular velocity of ip c (6t) is extremely small and the in- 
terference is not determined with de Broglie phase but 
this small phase. This leads an unusual neutrino inter- 
ference of the present work. The next singular term is 
from 1 1 A in & n j, and becomes Jg{x) / ^7ro- v \p v \ 2 , which is 
much smaller than Jgy) in the present parameter region. 
The magnitude is inversely proportional to \5t\ and is 
independent from the in 2 for the general form of wave 
packet also. 

Integrating t \ and ?2 in a finite T = T v - T„, we have a 
slowly decreasing term g(T, a> v ) and a normal term Go- 
g (T, o)y) is generated from the light-cone singularity and 
related term. g(T, oj v ) = g(T, cj v ) - g(oo, lj v ) where 



Tg(T, ajy) = -i 



I 

Jo 



\8t\ 



(15) 



and g(oo, co v ) = 0. The normal term, TGo, is from the 
rest. Due to the rapid oscillation in St, Go gets contri- 
bution from the microscopic 6t region and is constant 
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Angle between and P v [degree] 

Figure 1 : The angle between p n and p v dependence of the neutrino 
energy is given. The horizontal axis shows the angle and the vertical 
axis shows the neutrino energy at E n = 10[GeV]. The normal term 
has a value along the boundary and the diffraction term has a value in 
a broad area below the normal term. 



in a macroscopic T. This term does not depend on <j v 
and agrees with the normal probability obtained with 
the standard method of using plane waves. In the region 
Ipn-pv > fh, A„j(5x) does not have the light-cone sin- 
gularity and the diffraction term exists only in the kine- 
matical region 2p„-p Y < in. 

We compute a total probability next. From an inte- 
gration of neutrino's coordinates X v , the total volume V 
is obtained and cancelled with the normalization of the 
initial pion state. The total probability, then, becomes 
a sum of the normal term Go and the diffraction term 
g(T, w v ), 



d 3 p v Pn-p v {m 2 n - 2p„-p v ) 
(2tt) 3 E v 



[f(T, Wv ) + G ], 
(16) 



where A^3 = 8Tg 2 cr v and L = cT is the length of a decay 
region. At T — > oo, the diffraction term vanishes and 
the probability P agrees with the value of the standard 
calculation of plane waves. At a finite T, the probability 
has the diffraction component, g (T, <y v ), which is stable 
under variation of the pion's energy. 

Now we study each term of Eq. (1161 . In the normal 
term, Go, the energy and momentum are conserved well, 
and Go has a sharp peak at p„-p v = in 1 12. Hence the 
factor m 2 . - 2p n ■ p v in Eq. (TToT i becomes m 2 and a cosine 
of angle between p n and p v is determined uniquely by 
the relation, 



(ml 



■ m\ + 2E„E V ) 



2\p„\\p v 



(17) 



and the probability is proportional to mf. Integration of 
the neutrino's angle leads this integral independent of 
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the angle width, as far as it include the narrow peak. 
The value is independent also of cr v , which is consistent 
with the condition for the stationary state 11511 . 

The diffraction component, g(T, u> v ), on the other 
hand, is present in a kinematical region, \p v \{E n -\p n \) < 
Pn'Pv ^ m 2 12. The energy and momentum are not ex- 
actly conserved in the space-time dependent amplitude 
Eq. (TTOl l. Hence m 2 — 2p„-p v in Eq. ( TToT l is larger than 
m?, and the cosine of angle between p„ and p v is not 
uniquely determined. The dependence of the diffrac- 
tion term upon the angle between p„ and p v is presented 
in Fig.[T]for E„ - 10 [GeV]. The angle, which is deter- 
mined uniquely from the neutrino energy for the small 
wave packet, is not unique and is widely spread. From 
this behavior of the diffraction component, it might have 
been possible to reject the diffraction component and to 
take only the normal component in narrow band beam 
experiments if the neutrino's interaction position were 
observed. Comparison of the total cross section of sat- 
isfying this constraint with that of non-constraint events 
might have given direct signals of the diffraction com- 
ponent. 

The diffraction term is slowly varying with both the 
distance and energy. The typical length Lo of this 
universal behavior is Lo [m] = 2E v hclm 2 = 400 x 
£ v [GeV]/m 2 [eV 2 /c 4 ]. Its magnitude was found to be 
about 10-20 % of the normal term in a suitable exper- 
imental situation and depends on geometry of the de- 
tector. In this paper we analyze total neutrino nucleon 
cross section, which is connected with the neutrino flux 
and its high energy behavior is understood well from a 
standard quark parton model. 

5. Charged leptons 

A probability of observing a neutrino or a charged 
lepton in the pion decay at a position of a macroscopic 
distance is written in the form 



P = Pn 



In Eq. ( TT8l >, P„ onm i is the normal term that is obtained 
from the decay probability Go in Eq. (TToT l and P l di ff 
is the diffraction term that is determined from g in 
Eq. (TToT l. The latter has not been included in calculating 
the total cross sections before and its effect is estimated 
here. 

The diffraction term at a time T is described by its 
mass and energy in the universal form 



P'diff = Cdif/TgfJ, cod, 



(19) 




Figure 2: Values of g(T, oj) for v, e and /i at L = 10 [m], E = 1 [GeV] 
and m v c 2 = 1 [eV], m e c 2 = 0.5 [MeV] and m„c 2 = 100 [MeV]. 



is positive definite and decreases slowly with a distance 
L = cT, and vanishes at the infinite distance with the 
length scale 2cEi/m 2 . This scale is macroscopic for 
neutrinos but is 10~ 10 [m] or less for the electron and 
muon. The magnitude of g (T, a>i) is given in Fig. (0. 
At L = 10 [m], E = 1 [GeV] for the mass 1 [eV/c 2 ] (v), 
0.5 [MeV/c 2 ] (e) and 100 [MeV/c 2 ] (pi), the values are, 



g(T,(o Y ) ~ 3, 

&7>e)*0, 

0. 



(20) 



Hence the diffraction component at a macroscopic dis- 
tance is finite in neutrinos and vanishes in others. It 
is striking that the neutrino flux has an additional term 
and is not equivalent to that of the charged lepton even 
though they are produced in the same decay process. 
This diffraction term is generated by the tiny neutrino 
mass and the interference term generated with the su- 
perposed waves of forming the light-cone singularity. 
Because Pdtff is the interference term, it has unusual 
properties different from those of P„ orma i in the neu- 
trino flavour and the energy and momentum. The neu- 
trino diffraction furthermore is sensitive to the absolute 
neutrino mass. We study implications of the diffraction 
term to total cross sections, next. 

6. Total cross sections of v^-N scattering 

A total cross section of a neutrino nucleon scattering 
is written in the form, 

„ M N E v Gl 



■((2 + (2/3), 



(21) 



where C^iff is a constant and is obtained later. a>i are 
small in neutrinos and large in charged leptons. g(T, ai/) 



using integrals of quark-parton distribution functions 

rl rl 

q(x) and q(x), Q — J Q dxxq(x), Q — J Q dxxq(x). 
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The cross sectin is proportional to the neutrino en- 
ergy and a current value is cr v /E = (0.677 ± 0.014) x 
10~ 38 [cm 2 /GeV] |3[]. So experiments seem consistent 
with Eq. (f2Tb . However recent experiments of NOMAD 
J3l and MINOS ill] gave the total cross sections in wide 
energy ranges with small uncertainties and showed that 
the cross sections have slight energy dependences. They 
are compared with our theoretical calculations in the 
following. 

The diffraction term was identified only recently. 
Here we include the diffraction term into the neutrino 
flux. The neutrino flux in pion decays is given as a sum 
of the normal and diffraction terms in the form 



/ - fnormali^ + r diff)< 



(22) 



where rdiff is a ratio of the diffraction component 
over the normal component, and is a function of £ — 

m^L/lcEy, 



fdiff = d g(O, 



(23) 



where the coefficient do is determined from geometry. 

When the detector is located at an end of a decay vol- 
ume, the correction factor Eq. ( 1231 is used. In an actual 
experiment, the detector is located in a distant region 
from a decay volume. There is soil between them and 
pions are stopped in a beam dump. The neutrino diffrac- 
tion does not occur and neutrino propagates freely in 
this region. Since the wave packets of one cr v form a 
complete set, the wave packet size at the decay volume 
is the cr v determined with the detector. The neutrino 
flux at the end of the decay volume is computed with 
the diffraction term of the decay volume 's length L and 
the wave packet size of the detector. Wave packets of 
this cr v propagate freely from the end of a decay volume 
to the detector. The final value of the neutrino flux at 
the detector is found using the factor Eq. (l23l . When a 
neutrino changes flavour in this period, the final proba- 
bility for each flavour is written with a standard formula 
of flavour oscillation. 

A true value of the total cross section, a-{E) true , ob- 
tained with the total flux is connected with a cross sec- 
tion, cr(E) exp , obtained with only the normal component 
of the flux by the ratio 

o-(£)' rac = o-(E)" p — — . 

1 + r di f f 

Conversely the experimental cross section is written as 

a-(E) ex P/E = (1 + rw/KoiEf/E). (24) 

o-(E)"' ue /E is believed a constant so the E-dependence 
of ai_E) exp [E is due to the E-dependence of rdiff, 
Eq. (f23]l. 
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Figure 3: Total cross section of MINOS is compared with the sum of 
normal and diffraction terms. m v c 2 = O.leV 
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Figure 4: Total cross section of NOMAD is compared with the sum 
of normal and diffraction terms. m v c 2 = 0.2eV 

Since the diffraction term has different properties 
from those of the normal term, the corrections rdiff de- 
pends also on the geometry of the experiment and the 
material of the detector. In this paper we compute the 
cross sections using the experimental conditions of MI- 
NOS and MOMAD, which have presented precise cross 
sections of various energy ranges under the same condi- 
tion, and compare them with the experiments. 

The geometry of MINOS and NOMAD are the fol- 
lowing. The lengths between the pion source and the 
neutrino detector, L de t-so, and those of the decay re- 
gions, l^dec-reg, are: 

NOMAD : L det . so = 835 [m], L dec . reg = 290 [m], 
MINOS : Uet-so = 1040[m], L rfec _ reg = 675 [m]. 

The detector size is 7 x 3 x 3 [m 3 ] for both experiments. 

The wave packet size is estimated using the size of 
target nucleus. From the size of the nucleus of the mass 
number A, we have <x v = A 3 jm\. For various material 
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the value are 

o~ v = 5.2/m 2 ; l2 C nucleus, 

o~ v = 14.3/m 2 ; 54 Fe nucleus. (25) 

The total cross sections of MINOS and NOMAD ex- 
periments are compared with the theoretical values in 
Figs. (f3) and ©. In theoretical calculations, effects of 
a pion beam spreading is included by taking an average 
of initial pion's angle from to 10 [mrad] for NOMAD 
and from to 15 [mrad] for MINOS. An effect due to 
a finite size of the initial pion was also estimated. We 
found that an plane wave approximation which was em- 
ployed in this paper was very good. 

The cross sections decrease quite slowly with the 
neutrino energy, which may be difficult to understand 
with the standard theory. The theoretical cross sections 
obtained by including the diffraction component into the 
neutrino flux showed the same behavior and agreed well 
with the experiments. Two experiments are actually dif- 
ferent in the neutrino energy and geometry, but agreed 
with the theory. So the large cross sections at low energy 
regions may be attributed to the diffraction component. 

We have compared only NOMAD and MINOS here. 
Many experiments are listed in particle data 13J and 
most of them have similar energy dependences and 
agree qualitatively with the diffraction's presence. It 
is important to notice that the magnitude of diffraction 
component is sensitive to geometry and if a kinemati- 
cal constraint Eq. (fTTT i on the angle between p n and p v 
was required, only the events of the normal term was 
selected. Then the cross section should agree with that 
of the normal term. 

7. Summary and implications 

We showed that due to the diffraction component, 
Pdiff, the neutrino flux was modified. The total neutrino 
cross sections of NOMAD and MINOS agreed with 
the theoretical calculations obtained using the modified 
neutrino flux. Thus the existence of the neutrino diffrac- 
tion is consistent with experimental observations of the 
total cross sections. Although the diffraction component 
is determined by the average neutrino mass, its sensitiv- 
ity of the existing data is not sufficient to find an infor- 
mation on the neutrino mass below 0.3 [eV/c 2 ]. 

Other channels of interests are quasi-elastic or one 
pion production processes. The cross sections for v + 
n — > /j~ + p + (+7T°), v + p + — > yu~ + p + + n + and 
v + p + — > p + + n (+7T°) and the neutral current pro- 
cess are known well theoretically using models such as 



CVC, PCAC, and vector dominances and others. Re- 
cent experiments 11911 showed that the cross sections 
have excesses of 20 - 50% and are consistent with the 
diffraction terms. Especially a proton is known to have 
large wave packet size in matter due to its small mass 
and the diffraction due to a proton target is enhanced 
and gives a finite contribution despite its small energy 
ratio, m proton /^nucleus- 
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